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Abstract 

Complete chaotic synchronization of end lasers has been observed in a line of mutually coupled, 
time-delayed system of three lasers, with no direct communication between the end lasers. The 
present paper uses ideas from generalized synchronization to explain the complete synchronization 
in the presence of long coupling delays, applied to a model of mutually coupled semiconductor 
lasers in a line. These ideas significantly simplify the analysis by casting the stability in terms 
of the local dynamics of each laser. The variational equations near the synchronization manifold 
are analyzed, and used to derive the synchronization condition that is a function of parameters. 
The results explain and predict the dependence of synchronization on various parameters, such as 
time-delays, strength of coupling and dissipation. The ideas can be applied to understand complete 
synchronization in other chaotic systems with coupling delays and no direct communication between 
synchronized sub-systems. 

PACS numbers: 05.45.Xt, 05.45.-a, 42.65. Sf, 05.45.Vx 
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I. INTRODUCTION 



Synchronized chaotic^scillations have been found in many nonhnear systems, from lasers 
[2! to neural networks [3]. Many types of synchronization have been observed, including 
complete synchronization, phase-locking, and generalized synchronization, in the case of 
uni-directionally coupled systems. For an overview of the large body of work done on syn- 
chronization, and its sub-classes see for example 0, 0] . While extensive work has been done 
on mutually coupled systems, general analytic techniques for analyzing chaotic synchro- 
nization in time-delay mutually coupled systems have not been well developed. Essentially, 
there exist two methods for analyzing synchronization in coupled systems: A Lyapunov 
function approach jo] and a master stability approach Time delays considerably compli- 
cate the analysis, possibly introducing infinite degrees of freedom, and resulting in new types 
of dynamics [6;]. The present paper proposes an approach for understanding and predict- 
ing chaotic synchronization of time-delayed mutually coupled systems, possessing internal 
symmetry. 

In the case of internal symmetry in a system, where the equations of motion are invariant 
with respect to interchange of some variables, there is a solution where these variables are 
exactly equal [sl]. For example, for the case of three nonlinear oscillators coupled in a line, if 
the equations of the outside oscillators are identical, then given the same initial conditions, 
these oscillators will have identical dynamics, including the possibility of chaotic solutions. 
In that case, the system can be reduced to two coupled oscillators. If the symmetric solution 
is asymptotically stable to perturbations off the synchronization manifold, then the dynamics 
of outside oscillators are synchronized. Thus the requirement for synchronization in the long 
time limit is that the largest Lyapunov exponent, with respect to perturbations transverse to 
the synchronization manifold, is negative, resulting in decay of the initial perturbation back 
to synchronized state [7] . In general, Lyapunov exponents have to be calculated numerically. 
However, as will be shown, an analytic estimate can be made in some cases by linearizing 
about the synchronous state. 

Although there is extensive work on synchronization of coupled systems, studies of chaotic 
synchronization in time-delayed systems is much less extensive by comparison. Some analysis 
exists on synchronization of coupled semiconductor lasers without delays j^, Q]. However, it 
remains to be explained, for example, why in a coupled three laser system, outer lasers show 
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complete synchronization, in the presence of long time delays (compared to the internal 
dynamics of each laser) and no communication except via the middle laser, which itself is 
not synchronized with the end lasers (See Figure [2] below for an example). The present 
paper aims to explain this phenomena observed in lasers and other time-delayed systems 
using ideas from generalized synchronization. The paper is organized as follows: In Section 
II, the general equations for the three subsystems coupled by nearest neighbor interactions 
with delays are introduced, and the equations are linearized close to the synchronization 
manifold, using internal symmetry. Section III uses synchronization ideas developed in 
Section II to explain complete synchronization of semiconductor lasers in the presence of 
long delays and with no direct communication between the outer lasers. Section IV concludes 
and summarizes. 

II. LONG DELAYS AND GENERALIZED SYNCHRONIZATION 

We deal with mutually coupled, oscillatory time-delay systems possessing internal sym- 
metry, with respect to interchange of some variables. The system can be broken up into 
three coupled parts as depicted in Fig. [1) the "center" and two identical sub-parts, that 
possess symmetry with respect to interchange of variables: 



where Zi are vector variables of some dimension M and A^, for the outer and the middle 
subsystem, respectively. Variables zi and Z3 are symmetric with respect to interchange 
of variables, while the center system, Z2 may have different internal dynamics, given by 
F, and different coupling function G. The delay in the coupling terms is fixed and given 
by r, and the strength of coupling from the center to the outer identical subsystems by 
5i, while from the outer to the center by 82- Due to internal symmetry of the system, 
there exists an identical solution for the outside subsystems: zi(t) = Z3(t) = (j){t). If 
the conditional Lyapunov exponents calculated with respect to perturbation out of the 
synchronization manifold, Zi(t) = Z3(t) are all negative, then the two outer subsystems are 
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synchronized. This type of behavior where the two systems show identical dynamics, even 
chaotic ones, is called complete synchronization. Figure [2] shows an example of complete 
synchronization in the case of a mutually coupled three laser system using the model in 
Section III [10]. Calculating Lyapunov exponents is in general complicated due to the 
presence of time-delays in the equations. The coupling term containing delays, however, 
drops out if Eqs. (1), (3) are linearized about the synchronous state: zi(t) = zsit). To study 
the stability of the symmetric solution, we introduce new variables: Azi(t) = zi(t) — (j){t) 
and Az3(t) = zs{t) —(j){t). Linearizing transverse to the synchronization manifold, we have: 

ciAzi(t) 



dt 

dAzsit) 



J-Azi(t) (4) 
J ■ Az3(t) (5) 



dt 

where J is the M x M Jacobian matrix of partial derivatives evaluated at (f){t) 

J gF(0(t)) ^ 
dz 



(6) 



Here is the synchronous state that is determined by the dynamics of Eqs. (l)-(3) 
and the initial conditions defined on [— r, 0). Although the time delays and dependence 
on Z2 drop out of Eqs. (jlj) and ([5]), they are involved implicitly in determining (f){t), the 
synchronization manifold. Equations (jlj) and ([5]) are M-dimensional and therefore have M 
transverse Lyapunov exponents. The largest of them determines the stability of the trans- 
verse perturbation. So that the synchronized state zi(t) = Z3(t) = (j){t) is asymptotically 
stable, if Azi 3(t) ^ as t ^ oo or if all of the Lyapunov exponents in the linearized 
equations are negative. 

For 62 = 0, in Eq. (2), the dynamics of Zi 3 become that of a driven system, with Z2 acting 
as the driver. Then, the synchronized dynamics correspond to generalized synchronization 



ll| whereby the driven subsystem becomes a function of the driver: (j){t) = f(^*). Here, $ 



are the dynamics of the driver obtained by integrating Eq. (2), with 62 = 0: d^{t)/dt = 

While the exact form of the function between the driver and the driven systems can be 
rather complicated and difficult to obtain, its existence can be inferred from the synchroniza- 
tion of identical systems when started from different initial conditions and being exposed to 
the same drive. This method of detecting generalized synchronization using identical driven 
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systems is known as the auxiliary systems approach 12(]. In order for the driven subsystems, 
zi,3 to become synchronized, their dependence on initial conditions has to "wash out" as a 
function of time. This is due to the fact that dependence on initial conditions prevents syn- 
chronization by making the dynamics of Zi different from the dynamics of Z3. This "washing 
out" of initial conditions is provided by the dissipation in the system, which must therefore 
be either present in the coupling term, or in the uncoupled dynamics of the system itself. 
This can be seen by taking the sum of Lyapunov exponents, for Eqs. (jlj) and ([5]), which is 
related to contraction or expansion of the phase-space volume of the dynamics transverse to 
the synchronization manifold [3], 



M 

V A, = lim - In I det (*(Azi,3) (t)) |, (7) 

^ — ^ i^oo t 

i=i 

where ^ is the fundamental matrix solution to Equations (jl]) and (j5]). From Eq. ([7j), 
negative Lyapunov exponents sum corresponds to a contraction of phase-space of Azi 3 
dynamics as a function of time. In the case of synchronization in driven identical systems, 
all Lyapunov exponents transverse to the synchronous solution must be negative. Thus 
there is a contraction of phase-space to a single trajectory that is a function of the dynamics 
of the driver, This shrinking of phase-space is either caused by dissipative coupling or 
dissipation in the driven systems themselves. The effect of dissipation on synchronization 
can be illustrated using the case of a simple driven system, 

^ = -ea; + <l>(t) ^ = -ey + (8) 

taking the difference between the two variables, A = x — we get, d/\/dt = — eA, where 
t) in Eq. ([8]) is the chaotic signal provided by the driver (it can also be a noisy signal) 



14l |. Thus the difference in initial conditions between x and y decreases at the rate of 



dissipation, e, leading to synchronization for large times. 

In mutually coupled systems, 62 7^ 0, the dynamics of Z2 are affected by zi and Z3. In this 
case, the synchronized state, 4>{t), may depend on the initial conditions of all of the three 
subsystems, {zi,Z2,Z3}, so that (f){t) can not be the result of generalized synchronization, 
in a strict sense. However, it takes a time interval of 2r for any change in the dynamics of 
systems zi^s to affect the motion of these systems via mutual coupling. During this time 
interval of length, 2r, zi 3 can be viewed as driven by z^, since the signal zi 3 receives during 
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that time-interval is not affected by its dynamics on that intervaL Therefore, we examine 
the dynamics in a time period on the order of the delay time , r. 

The initial dynamics of zi s{t) on the time interval < t < t^ + T, affect the dynamics of 
Zi,3(t) on the time interval to + 2T < t < Iq + St via the mutual coupling term, G (z2 (t — r)) 
(see Eqs. (1) and (3)). For chaotic systems, the trajectories are not correlated in time, and 
we assume 



for i = {1, 2, 3} and to significantly longer than the average period of oscillation. The above 
equation is true, in general, for non-periodic oscillations. This can be seen by expanding 
the signal in a Fourier series: Zi{t) = Yl'^=-oo ^nCos [nt + (p). Then, if there is a significant 
periodic component in Zi{t) of amplitude An, Eq. ([9]) will be proportional to Af^ whenever 
to is a multiple of 27r/n. Thus, for example, Eq. ([9]) may not hold if the amplitude of an 
optoelectronicly coupled laser (to be discussed in the following section) is too close to the 
threshold, where the behavior can be approximated as coupled linear oscillators (see Eqs. 
f fT0|) - f|T2|) ). leading to a significant regular oscillatory component in the signal. In this case, 
there may be a resonant interaction between the lasers, which is sensitive to the specific 
value of the coupling delay, r. For chaotic dynamics, we can substitute to = 2r into Eq. (Q, 
where r is the delay, to see that there is no significant correlation between the dynamics on 
that time scale, so that over the round-trip time of 2r, the identical sub-systems zi 3 can be 
viewed as driven by some uncorrelated chaotic signal coming from Z2. This assumption of a 
driver is only strictly valid on the time interval within the round-trip time, since for longer 
time intervals the initial conditions of the outer lasers at the beginning of the interval will 
affect their dynamics, via the middle system, at a later time within the interval. 

By perturbing the dynamics of the outer subsystems from the synchronized state, it can 
be shown that complete synchronization of the end subsystems in the presence of long delays 
is similar to generalized synchronization, where the middle subsystem acts as the driver for 
the outer ones. After the symmetric subsystems synchronize, zi_3(t) = 0(t), one of them 
can be suddenly perturbed from its symmetric state to an arbitrary position in phase-space 
at some t = ti. In that case, the perturbed system, as well as the unpertubed one, will 
receive the exact same signal from Z2 as before, for t < ti + 2r. If the systems synchronize 
again at some point during t < ti + 2r, we will again have zi 3(t) = 0(t), where 0(t) has 
not been affected by the perturbation during that time interval. Thus the synchronized 




(9) 
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state, (f>(t) = zi^sit) is clearly independent of perturbations of subsystems zi^a and must 
therefore be some function of the middle subsystem, Z2. This however is the same as what 
happens in generalized synchronization, with the difference that the trajectory of itself 
may be affected by the initial starting conditions of the symmetric subsystems. The following 
section uses a system of semiconductor lasers as an example for application of these ideas 
to understand complete synchronization of the end lasers in a three laser system. 



III. SYNCHRONIZATION OF SEMICONDUCTOR LASERS WITH DELAYS 



The discussion of the previous section can be applied to the study of synchronization of 
a three laser system with delays. A schematic diagram is shown in Fig. ([1]), where the outer 
lasers are identical, while the middle laser is detuned. 

The scaled equations of coupled semiconductor lasers have the following form 



15 



16|: 



— = xi [l + yi) 
— = -yi - exi (ai + hiyi) + 52?/2 {t - r) (10) 

^ = /? (-1/2 - e/3x2 (as + 62Z/2)) + 5i (Vi (t - t) + {t - r)) (11) 

dy3 . X 

^ = X3(l+y3) 

^ = -Z/3 - ea;3 (ai + hyi) + 52y2 it - r) (12) 

Eqs. fllOI ) -( fT2l) have the same form as Eqs. (l)-(3) with Zi = {yi, Xi}, where 2 = 1,3 and i = 2 
for the outer and middle lasers, respectively. Variables y^ and Xi denote scaled intensity and 



inversion of the ith laser, {ai, 02? ^i? ^2} are loss terms, and e is the dissipation. (See 17| for 
details of the derivation form the original physical model.) The above equations are coupled 
via laser intensities, y^, using optoelectronic incoherent coupling that does not contain phase 
information, unlike the coherent coupling in Fischer et al pLj]. Previously, the dynamics 



of two electronically coupled lasers have been explored 18|, showing lag synchronization 
for the case of two lasers, and isochronal synchronization if feedback is added [3]. The 
above equations are scaled so that the relaxation frequency is equal to one. In the typical 
experimental set-up, the relaxation oscillations are on the order of 2 — 3 ns. Since the delay 



time is scaled by the relaxation frequency, and is at least an order of magnitude higher 
(for long delays), a typical delay time used in simulations could be about r = 60, which 
corresponds to about 20 — 30 ns. 

In the absence of dissipation, the uncoupled system, 5i = 82 = 0, is a nonlinear conser- 
vative system, with behavior similar to a simple harmonic oscillator for small amplitudes. 



and becoming more pulse-like at high amplitudes [17!]. Dissipation, however, leads to energy 
loss, so that in the absence of coupling between lasers, the dynamics would settle into a 
steady state. Thus mutual coupling acts like a drive by pumping energy into the system. 
For most cases, it can be assumed that dissipation is small: e <^ 1. Detuning of the middle 
laser from the outer ones is given by (3. 

The system described by Eqs. f lTU]) - f[T^ shows complete synchronization of outer lasers 
over a whole range of parameters. Figure [2] shows that while the outer lasers can become 
completely synchronized, there may be no apparent correlation between the middle and the 
outer lasers. Since the outer lasers are identical, there is a solution of Eqs. (fTOl) - ( fT2l) where 
2/1 = 2/3 = Y{t) and xi = xz = X{t). In this case, Eqs. (ITOl) - (|T2l) reduce to four differential 
equations. The solution yi = and Xi = X3 is stable if the Lyapunov exponents transverse 
to the synchronization manifold are negative. To investigate the stability of the synchronized 
state we linearize about the synchronous solution (pit) = {X{t),Y(t)}. Applying Eqs. (j4])- 
([6]) to Eqs. ([10]) and ([12]), we get 



(13) 



Axi,3\ _ (-e[a^ + b^Y{t)] -[1 + eb^X {t)]\ / Ax(t)i,3 
Ai/i,3y \ l + Y{t) X{t) ) \^Ay(t)i,3 

where {Ax(t)i,3, Ay(t)i^3} are perturbations of outer oscillators from the synchronous state 
{X{t),Y{t)}. This synchronous state is obtained by starting the outer oscillators from the 
same initial conditions and perturbing the system at some time, t = ti. The perturbation 
will not affect the coefficient matrix in Eq. f[T3l) until t > ti + 2r. So that in the time 
interval of 2r the dynamics off the synchronization manifold can be viewed as driven by 
an uncorrelated chaotic signal {X{t),Y{t)}. We can now apply Abel's formula to Eq. 
f[T3]) . which relates the Wronskian of the linearized system to the trace of the matrix [20 1. 
Dropping the subscripts on linearized variables, we get. 



W (t) = det 



Ax Ay 
Ax Ay 



exp (^j^ {X{s) - e (ai + biY{s))} ■ ds^ (14) 
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The Wronskian gives the phase-space volume dynamics of the system {Ax(t), Ay(t)}. Equa- 
tion f|T^ is vahd over the integration interval of twice the delay: ti < t < ti + 2r. This 
is due to the fact that it takes a time interval of 2r for a perturbation in the outer laser 
to affect its dynamics via mutual coupling from the middle laser. Thus, during the time 
interval of 2r, the perturbed system acts like a driven system in that its dynamics do not 
affect the signal it receives, and therefore do not change the synchronized state, {X,Y}, 
making it independent of {Ax, Ay} dynamics over the integration interval. 

Since the variable Y{t) is the scaled intensity of the laser, from Eqs. (fT0l) - (fT2l) . its 
minimum possible value is —1. Thus for ai > bi (a typical case), the contribution of the 
dissipation term to the Wronskian is always negative. The variable X{t), on the other hand, 
is symmetric about zero, and thus averages out to zero when integrated over a single period 
of oscillation. It follows that if the integral in Eq. f|T^ is taken just over a single oscillation 
of the laser, we get 

/ {X{s)-e{ai + biY{s))}-ds = -e{ai + biY)T <0 (15) 
Jti 

where T is the period of a single oscillation, and Y is the average of Y over a single period 
(unlike X, the Y variable is not symmetric about zero, which can readily be seen in the 
pulse-like fluctuations of lasers at high intensities). It follows, that X{s) in Eq. f|T^ averages 
out to zero if the integral is done over many periods of oscillation, while the dissipation term, 
multiplied by e, provides a continuous negative component. If that continuous negative com- 
ponent builds up sufficiently over the integration interval to overcome any fluctuations in 
X{s), we then have a continuous shrinking of the phase-space of perturbed dynamics, indi- 
cating synchronization. Integrating the exponential term in Eq. f|T4|) over many oscillations 
and using Eq. (fT5|) . we get, 

[ {X{s)-e{ai + biY{s))}-ds^-e{ai + biY){t-ti)+ [ X{s)-ds (16) 

Jti Jti+nT 

where t — ti is the total integration interval, Y is the average value of intensity over that 
interval, and n on the integration limits is the total integer amount of full oscillations that 
fit into the integration period: t — ti — nT < T. Here, T is the average period of oscillation 
over the integration interval. Thus the integral of X{s) on the right hand side is only over 
a single uncompleted oscillation. This integral, however, may still be significant compared 
to the e term, since its fluctuations are comparable to X{s), because the integration period, 
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T, is of order unity (due to scaling in the equations), while e multiplying the other term 
is small. It follows that sufficiently long integration times, t — ti, are required in order for 
the dissipation term to dominate. We can now set an upper bound for the integral on the 
right-hand side of Equation fITB]) . 

/ X{s)-ds<7i\X{t)U. (17) 

where |X(t)|m(i2. is the maximum fluctuation of inversion over the interval of twice the delay 
time. The above bound may not be valid at energies far above the threshold, when the 
laser behavior becomes pulse-like with a period that is significantly longer than the scaled 
relaxation period of 27c. This may be another reason why there is a loss of synchronization 
at higher coupling strengths, which lead to higher amplitudes of oscillation, with lower 
frequencies. Requiring Eq. (fT6l) to be less than zero and using Eq. (ITTI) . we can now 
obtain a bound above which the dynamics tend toward the synchronization manifold over 
the interval of twice the delay time, 

-Te(ai + biY) > \X{t)\^ax (18) 

TT 

where we have used 2t = t — ti for the integration interval. The above inequality insures 
the right-hand side of Eq. ( |T6l) is negative over the interval of twice the delay time. This 
in turn ensures the shrinking volume of transverse phase-space dynamics given by Equation 

m- 

Equation f|T8l) gives a condition for the phase-space volume of transverse dynamics 
to contract over the interval of twice the delay. For sufficiently long delays, where 
re (ai + feiF) ^ \X(t)\max, Eq- (fT4l) can be approximated as 

In {W{t)) = In \AxAy - AyAx\ ~ - /" e (ai + &iF(s)) ■ ds. (19) 

Jti 

where a natural log of W{t) was taken. The above equation is a monotonically decreasing 
function of t. This means that the phase-space volume of the system perturbed from the 
synchronization manifold contracts as a function of time. Applying Eq. ([7j) to Eq. f|T9l) . 
where det (^(Azi 3)) = W{t), the sum of transverse Lyapunov exponents for the dynamics 
off the synchronization manifold described by Eq. ( |T3l) can now be approximated as, 

Ai + A2 ^ -e (ai + h?) (20) 
10 



Equation (120|) indicates that for sufficiently long integration times (requiring sufficiently long 
delays), the sum of the Lyapunov exponents should be negative, indicating the shrinking of 
phase-space volume of dynamics transverse to the synchronization manifold. 

Figure [3] shows numerically computed sum of Lyapunov exponents, and corresponding 
correlations of the outer lasers, as a function of dissipation, e, for two values of the delay, 
T = 120 and r = 240. The fluctuations in the sum of Lyapunov exponents correspond well 
to the fluctuations in the correlation function of the outer lasers, with desynchronization 
when the Lyapunov sum increases above zero. As might be expected from Eq. (ITS!) . longer 
delays mean synchronization at lower values of dissipation, since the dissipation term in the 
exponential in Eq. (1141) dominates for sufficiently long delays. Increasing r by a factor of 
two, however, does not decrease the bifurcation value of e for the onset of synchronization 
by a factor of two, as might be expected from Eq. ffTSl) . This is probably due to the 
decrease in fluctuations, \X{t)\max, as the dissipation in the system is increased, leading to 
synchronization at a lower value of e than might otherwise be expected. After Eq. ( fTSl) is 
satisfied, resulting in the onset of synchronization, the sum of Lyapunov exponents has a 
negative linear dependence, given by Eq. ( l20l) . This is in agreement with Fig. [3] which shows 
this negative linear dependence of Lyapunov sum on dissipation, with a slope of around —2.6, 
a reasonable value for the parameters used of ai = 2, 6i = 1 and intensity, F ~ 1. In general, 
the average intensity of the dynamics, Y, depends on the coupling strengths, 6i, and 

While Eqs. ( IT9l) and ( |20l) predict the shrinking of phase space for the dynamics transverse 
to the synchronization manifold, to guarantee stability both Lyapunov exponents have to 
be negative, or the solution will blow up along the unstable direction. To find out whether 
the synchronous state is stable, consider again the Wronskian, W{t) = Ax Ay — Ay Ax. 
Substituting for Ax and Ay from Eq. f|T3l) . we get 

W{t) = {1 + Y{t)} ■ (Ax)' + {1 + ehX{t)} ■ {Ayf + {eK + hY{t)] + X{t)} ■ Ax Ay (21) 

For |e6iX(t)| < 1 (a reasonable assumption since e <^ 1), terms quadratic in Ax and Ay 
are always positive, indicating rotation. In Eq. (fT9l) . W{t) is a monotonically decreasing 
function of time, with W{t) ^ Q as t ^ oo. Therefore both Ax, Ay ^ as t ^ oo, due 
to the presence of positive terms quadratic in Ax and Ay in Eq. ( 1211) . It follows that for 
sufficiently long delays in the system, the synchronized state is stable, and therefore all the 
Lyapunov exponents transverse to the synchronization manifold are negative. The stability 
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of synchronized state is due to cross-terms in the matrix in Eq. (fT3|) . which come from 
rotation, leading to the nonhnear exchange of energy between inversion and intensity of the 
laser. This rotation introduces positive quadratic terms in Ax and Ay into Eq. (^T\i and 
leads to the spiraling of the phase-space volume towards zero, rather that blowing up along 
one direction, while shrinking along another. 

Figure shows the sum of Lyapunov exponents as a function of delay. The Lyapunov 
exponents are negative for all r > 170, (corresponding to about 60 ns) resulting in complete 
synchronization of the outer lasers, as shown in Fig. |3)d. At the same time, the outer lasers 
are not synchronized with the center one. Fig. Hb. The fluctuations in correlations of the 
outer lasers match well the fluctuations in the Lyapunov sum, with correlations increasing 
whenever the Lyapunov sum decreases. Figure H] agrees well with the above analysis, since 
sufficiently long delays (see Eq. f|T8l) ) are needed for the Lyapunov exponents to become 
negative, leading to synchronization. After the onset of synchronization, Eq. (l20l) becomes 
valid, so that the Lyapunov sum becomes independent of delays. This is confirmed by 
the straight horizontal line in the figure, after the outer lasers synchronize. The degree of 
synchronization is given by the correlations function. 

From Eq. fl20l) . the negative Lyapunov exponents, leading to stability of synchronous 
state, are the result of dissipation, e, in the end lasers. This is to be expected since mutual 
coupling pumps energy into the system, as can be seen in Eqs. ffTOj) - f|T2|) . Therefore, 
some dissipation in the outer subsystems themselves is essential in order to "wash out" their 
dependence on initial conditions and make them a function of the dynamics of the middle 
laser, as would be required in the case on complete synchronization. 

The amplitude of laser oscillations depends on the coupling strengths, 6i and 62, as well 



as the dissipation. It was shown 



151 ] that in the case of a two laser system there is a 



bifurcation value for the onset of oscillations that is a function of product of the coupling 
strengths, 6162- Increasing the coupling strengths increases the role of nonlinearities in the 
system and the intensity of laser oscillations. Thus for low values of the coupling strengths, 
the dynamics given by Eqs. (ITOl) - ( |T2l) are more regular. At low intensities, the dynamics 
of individual lasers are close to that of a simple harmonic oscillator, as can be verified by 
substituting low values of {xi,yi} into Eq. (fTOl) . for example. In order for Eq. ([9]) to be 
valid, the dynamics have to be uncorrelated over the time interval of the delay. Thus the 
equations derived in this section are valid for chaotic regime which requires sufficiently high 
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product of coupling strengths, 5i52- In this case, it can be assumed that the outer lasers 
are driven by an uncorrelated signal from the middle laser over the time interval of 2r. Any 
synchronization on that interval would then be analogous to generalized synchronization 
that occurs in a uni-directional system, with the middle laser acting as the driver for the 
outer ones. Since increased coupling pumps more energy into the system, thereby increasing 
the effect of nonlinearities, the Lyapunov exponents may increase above zero, leading to 
desynchronization of the outer lasers. In this case, longer delays in coupling may be required 
in order for the outer lasers to synchronize. This effect is illustrated in Fig. [5], which shows 
the sum of Lyapunov exponents as a function of coupling strengths for two different delays, 
r = 60 and r = 120. There is an abrupt increase in Lyapunov exponents above zero, due to 
increased nonlinearity, as the coupling strength is increased. Increasing the delay however 
to r = 120 leads to synchronization for a greater range of coupling strengths, as compared 
to r = 60. The corresponding correlations as a function of coupling strengths are shown 
in Fig. [6l Desynchronization at higher coupling strengths, and the synchronizing effect of 
increased delays is in agreement with Eq. (|T8l) . Since higher coupling strengths lead to 
greater fluctuations in X{t), longer integration times are required in order for Eqs. ( |T9i) and 
( 120|) to be valid, leading to synchronization at longer delays, r. 

Figures [3] - [6] show that Eq. (!20l) correctly predicts the independence of Lyapunov sum 
on delays and coupling strengths and a negative linear dependence on dissipation, once 
synchronization sets in. Synchronization, on the other hand, occurs once the condition 
expressed in Eq. ( IT8|) is satisfied, leading to the continuous shrinking of the phase-space 
dynamics transverse to the synchronization manifold. 

IV. CONCLUSION 

Ideas from generalized synchronization were used to explain complete chaotic synchro- 
nization of mutually coupled systems in the presence of long delays. Since identical outer 
subsystems synchronize due to a common input from the middle subsystem, complete syn- 
chronization is similar to the one occurring in the auxiliary system set-up, with the exception 
that all subsystems are mutually coupled. This leads to the dependence of common input 
to the outer subsystems on history of the dynamics. Complete chaotic synchronization is 
the result of the outer systems becoming a function of the middle one, as would happen in 
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the case of generalized synchronization. 

Due to the symmetry of the outer subsystems, the dynamics hnearized about the synchro- 
nization manifold are independent of explicit coupling. Transverse Lyapunov exponents can 
then be calculated to determine the stability of the synchronous state. Since over the time 
scale of twice the delay interval, the outer subsystems can be viewed as driven by a common 
chaotic signal from the middle subsystem, the analysis is considerably simplified, allowing 
for calculation of phase-space volume dynamics transverse to the synchronization manifold. 
The transverse phase-space volume dynamics were analyzed for the case of three mutually 
coupled semiconductor lasers. It was found that for sufficiently long delays, the synchronized 
state is stable. The sum of Lyapunov exponents transverse to the synchronization manifold 
was found analytically and shown to have a negative linear dependence on dissipation, in 
good agreement with numerical calculations. This also confirmed the intuition that synchro- 
nization is the result of dissipation, e in local dynamics of the lasers themselves, since the 
coupling between lasers is not dissipative. The analysis also explains the effect of various 
parameters on synchronization, such as coupling strengths, delay time, and dissipation, and 
is supported by numerical simulations over a range of parameter values. Thus, it was shown 
analytically and confirmed numerically that after the onset of synchronization, the stability 
of the synchronous state (as given by Lyapunov exponents) depends linearly on dissipation, 
but is independent of the delay time and coupling strength. 
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Laser 1 Laser 2 Laser 3 



FIG. 1: A schematic showing how three lasers are coupled in a line. The outer two lasers (circles) 
are identical, while the middle laser (square) is detuned from the rest. 
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A mutually coupled delay laser. 
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FIG. 2: Top: Intensity of Laser 1 vs. Laser 2. Bottom: Laser 1 vs. Laser 3. Straight line indicates 
complete synchronization of outer lasers, r = 30, e = VO.OOl, 5i = 82 = 6.5e 
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A mutually coupled delay laser-T=120. 




FIG. 3: a) Sum of Lyapunov exponents as a function of dissipation, e, for r = 120. b) Corresponding 
correlations between outer lasers, r = 120. c) Sum of Lyapunov exponents vs. e, for r = 240. d) 
Corresponding correlations between outer lasers, r = 240. In all cases, ai = a2 = 2, 6i = 62 = 1; 
61 = 82 = 0.2, P = 0.5. 
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FIG. 4: a) Numerically computed sum of Lyapunov exponents as a function of delay, r. b) 
Corresponding correlations of outer lasers, c) Correlations of the middle and outer lasers, shifted 
by the delay time to maximize correlations, e = \/0.001, 6i = 62 = 7.5e, (3 = 0.5. 
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FIG. 5: a) Sum of Lyapunov exponents as a function of coupling strength, 6i = 82, for r = 60. b) 
r = 120. e = VO.OOl, /? = 0.5. 
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FIG. 6: Correlations corresponding to Fig. [5j a) Correlation between the middle and one of the 
outer lasers, r = 60. b) Correlations of outer lasers, r = 60. c) Correlation between the middle and 
one of the outer lasers, r = 120. d) Correlations of outer lasers, r = 120. Outer lasers synchronize 
for greater range of coupling strength as the delay is increased. The middle and the outer lasers 
show little correlation for all values of the coupling strengths. 
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